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Abstract 

Calculations of spectral darkening and visibility functions for the brightness oscil- 
lations of the Sun due to global solar oscillations are presented. This has been done 
for a broad range of the visible and infrared continuum spectrum. The procedure for 
the calculations of these functions includes the numerical computation of depth- 
dependent derivatives of the opacity caused by p modes in the photosphere. A 
radiative transport code was used for this purpose in order to get the disturbances 
of the opacities due to temperature and density fluctuations. The visibility and 
darkening functions are obtained for adiabatic oscillations under the assumption 
that the temperature disturbances are proportional to the undisturbed temperature 
of the photosphere. The latter assumption is the only way to explore any opacity 
effects since the eigenf unctions of p mode oscillations have not been obtained so far. 
This investigation reveals that opacity effects have to be taken into account since 
they dominate the violet and infrared part of the spectrum. Due to this, the visibility 
functions are negative for those parts of the spectrum. Furthermore, the darkening 
functions show a wavelength dependent change of sign for some wavelengths due 
to these opacity effects. However, the visibility and darkening functions under the 
assumptions used are in contradiction with observations of global p mode oscilla- 
tions. But it is beyond any doubt that the opacity effects influence the brightness 
fluctuations of the Sun due to global oscillations. 

Keywords: Oscillations, Solar; Waves, Acoustic; Integrated Sun Observations; He- 
lioseismology, Theory; Spectrum, Continuum 



1. Introduction 

Since the ACRIM experiment (Woodard, 1984) it is well known that solar oscillations 
lead to brightness fluctuations of the Sun as a star. However, the spatial integration 
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of the flux disturbances over the whole solar disk suppresses high frequency modes 
and photometric observations of the Sun as a star can only resolve spatial oscillation 
modes with low harmonic degree (I < 3). But if the visibility of these low degree 
modes is suppressed at some positions on the solar disk due to opacity effects, then 
it might be possible to observe oscillation modes with harmonic degree I > 3. 

The space experiment IPHIR onboard the PHOBOS spacecraft (Frohlich et al., 
1988) showed that the amplitudes of the brightness fluctuations due to p modes 
are different for different optical wavelengths. Zhugzhda, Dzhalilov, and Staude 
(1993) showed with the consideration of non-adiabatic waves in a uniform, non- 
grey atmosphere that this problem cannot be solved with the introduction of the 
Rosseland mean opacity. With the calculation of visibility functions for low-degree 
non-radial oscillations and a comparison with IPHIR data, Toutain and Gouttebroze 

(1993) suggested two limiting cases for the explanation of the origin of the brightness 
fluctuations: 

• intensity and flux perturbations, caused by adiabatic oscillations, are driven 
by opacity effects and thus have to be taken into account for calculations of 
visibility functions 

• the visibility functions depend on non-adiabatic effects and opacity disturbances 
can be neglected (blackbody approximation). 

The direct way to solve the problem is to solve the eigenproblem for nonadiabatic os- 
cillations for the solar model which includes a standard model of the photosphere. To 
our knowledge nobody has done this so far since the set of integro-differential equa- 
tions of radiative hydrodynamics has to be solved. Staude, Dzhalilov, and Zhugzhda 

(1994) and Zhugzhda, Staude, and Bartling (1996) proposed to use another approach 
which assumes that the temperature and density fluctuations due to p modes are 
given. Under some special assumptions this makes it possible to obtain so-called dark- 
ening functions which shows the dependence of brightness fluctuations with respect to 
the position on the solar disk. The intention was to develop a qualitative explanation 
of brightness oscillations in order to improve observations of solar oscillations. The 
different wavelengths under consideration indicated that the brightness fluctuations 
of the most dominant oscillation mode (I = 0) vanishes for certain positions on the 
solar disk, e.g. the darkening function changes its sign. With the assumption that 
these calculations can be easily extended to modes with harmonic degree < £ < 3 
and that those functions will also show a change of sign, it is seen as an indication 
that modes with higher harmonic degree, £ > 3, might be observable. 

Zhugzhda, Staude, and Bartling (1996) carried out calculations of darkening func- 
tions for a few wavelengths only. Besides, the visibility functions were not calculated 
at all. In this paper the calculations of darkening and visibility functions are per- 
formed for a wide range of wavelengths, from the near UV to the infrared in order 
to see whether the opacity effects are essential for this range of the spectrum. This 
investigation will help to analyze the observations from photometric instruments like 
DIFOS on board CORONAS-F (Lebedev et al, 2004). 

2. Intensity and Flux Fluctuations due to P Modes 

To avoid the full solution of the equations of radiation hydrodynamics (Zhugzhda, 
Staude, and Bartling, 1996) we expand the general expression for the intensity in a 
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series around the equilibrium state up to the first order of small disturbances of the 
opacity. Assuming local thermodynamic equilibrium (LTE), the fluctuations of the 
emergent intensity (SI V ) normalized to the specific intensity (Iov) at the center of 
the solar disk are (Zhugzhda, Staude, and Bartling, 1996) 

Sly _ f°° e -T„/n 
hv Jo 



-1 

(1) 

where t v is the optical depth at the frequency v, B v the Planck function, p = cos 9, 
and 6 the polar angle. The unperturbed, outgoing intensity, originating at optical 
depth tv, is 

I v (t v ,») = e-^-^^B v {r' v )^. (2) 

The fluctuation of the opacity (5k v (p,T)) is given by the fluctuation of the density 
(Sp) and the temperature (ST), 

Sk^ dlnKo^ ST ^ dhiKoiy Sp 
kou dlnT To dlnp po 

Equation (1) considers the various contributions to the intensity fluctuations. The 
first term appears because of the temperature fluctuation in the source function, 
which can be described by the Planck function in the case of LTE. The second term 
appears because of the perturbation of the opacity due to disturbances in density and 
temperature. An approximation for the temperature and the density perturbations 
will be used. The solution yields the intensity fluctuation of a non-grey atmosphere 
caused by small disturbances. This can be used for the disturbances due to p modes 
since they are small enough and the linear approximation works perfectly. 

To simplify Equation (1) it is assumed that disturbances of the equilibrium state 
are adiabatic. In this special case the temperature and density disturbances are 
connected by a simple relation: 

f = (r 3 -i)^, (4) 

where is the third adiabatic exponent. For a neutral or fully ionized gas it is 
constant and has the value — 1 = 2/3, whereas for a partially ionized gas it can 
reach values of T3 — 1 < 0.1. The first case is given for the photosphere, which is 
completely dominated by neutral hydrogen. The second case applies for the deeper 
photosphere and the chromosphere where the gas is partially ionized. 
In the case of adiabatic disturbances, the intensity fluctuations depend only on 
temperature fluctuations and on parameters of the quiet undisturbed photosphere. 
The dependence of the intensity fluctuations on p is a darkening function for the 
adiabatic disturbances. To obtain the darkening function for p modes the depen- 
dence of the temperature disturbances on depth (ST(t v )) has to be substituted into 
Equation (1). But this function is not known, even for adiabatic disturbances. In 
order to make a qualitative analysis of the opacity effects we simplify the problem 



AB u (t v ) ST Sk v 

T] r 7fT H (B v (Tu) - h(T v ,p)) 

d In I o Jo nov 



d-rv 



e Tv B v (t v )&t v 



SolPhys_Nutto.tex; 8/03/2008; 9:14; p. 3 



C. Nutto et al- 



and it is assumed that the disturbances are proportional to their undisturbed values 
(Zhugzhda, Staude, and Bartling, 1996), 



5T 



constant 



(5) 



which is a rather crude approximation. It holds only for continuum radiation which 
originates from a thin layer of the solar atmosphere. With this assumption the relative 
disturbances 5T/Tq and Sp/po are independent of the optical depth (tv) and can be 
factored out of the integral. Now, with relationship (4) the relative intensity variation 
can be described in units of 5p/p or 5T/T. This leads to the darkening function 



A, 
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(6) 



which in our case is expressed in units of ST/To. 

To obtain the visibility function, the darkening function has to be integrated over 
the entire solar disk. This yields the following expression for the flux perturbation 
for the most dominant mode £ = 0: 
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3. Numerical Calculation of Darkening and Visibility Functions for 
Solar Oscillations with Degree £ = 

For the evaluation of the darkening functions one needs to know radiative transport 
quantities for the solar atmosphere as functions of the optical depth t v . For this 
purpose the radiative transport code RH by Han Uitenbroek was used (Uitenbroek, 
2001). The code is based on the MALI (Multi-level Approximate Lambda Itera- 
tion) formalism by Rybicki and Hummer (1991). It solves the combined equations 
of statistical equilibrium and radiative transport under the general assumption of 
non-LTE conditions for a multilevel atom in a given plasma. In our case we choose 
hydrogen to be the considered atom. In addition to opacities and emissivities from 
the transitions in the hydrogen atom, the code accounts for background radiation 
sinks and sources due to other atoms, molecules, and all relevant continuum processes 
such as H~ bound-bound and bound-free processes, Rayleigh scattering off neutral 
hydrogen and Thompson scattering off free electrons, helium, H2 and hydrogen free- 
free processes. The considered atoms and molecules for the background radiation and 
opacities are listed in Table 1. In addition, it is necessary to provide the radiative 



1 http: / /www. nso.edu/staff/ uitenbr/rh.html 
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Table 1. Atoms and molecules that are considered as background 
elements by the radiative transport code. 



Atoms He, C, N, O, Na, Mg, Al, Si, S, Ca, Fe 

Molecules H 2 , H+, C 2 , N 2 , 2 , CH, CO, CN, NH, NO, OH, H 2 



transport code with a model of the solar atmosphere. We choose the semi-empirical 
one-dimensional quiet-Sun model FAL-C by Fontenla et al, (1999). 

For the calculation of the darkening function we are not only interested in the 
radiative quantities alone, but also in their disturbances. Considering Equation (1), 
we need the disturbances of the Planck function due to temperature variations and 
the disturbances of the opacity due to variation of temperature and density. 
In order to calculate the derivatives 

dB v \ (d\riK v \ , fdhiK v \ 

and > ( 8 ) 



SlnT J ' V SlnT J ' V d\np 

we disturb the quantities T and p of the atmospheric model which is then used as 
input for the radiative transport code. We use temperature disturbances of —5 K 
and +5 K, which is less then one percent of the undisturbed temperature of the 
atmosphere and the derivatives behave linearly under these small temperature dis- 
turbances. With Equation (4) for the adiabatic relation, it can be assumed that the 
relative disturbances of the density are around 1%. Thus, the hydrogen and electron 
populations are modified accordingly in the model of the atmosphere. While one 
modification is applied, all other values are kept constant. This procedure allows us 
to calculate the partial derivatives of the Planck function and the opacity using the 
centered-difference approximation (Pozrikidis, 1998). 

Our interest is in the calculation of the darkening and visibility functions for 
a broad range of the visible and infrared spectrum. This has been done for 79 
wavelength points between 250 nm and 2100 nm. We have to avoid any lines since we 
assume a temperature disturbance that is constant, see Equation (5). This is only 
justified as long as the radiation originates from a thin layer. This applies for the 
continuum radiation, but spectral lines form in higher layers and are not constrained 
to a small part of the solar atmosphere. 



4. Results 
4.1. Derivatives 

First, the derivatives of the Planck function [dB v /dhxT) were calculated. The 
Planck function has the advantage that the derivative can be evaluated analytically 
and the numerical values can be compared with the analytical function. A comparison 
showed that the numerical calculation of the derivatives coincides with the analytical 
values within a few tenths of a percent. 

In Figure 1 the derivative of the opacity with respect to temperature is shown as a 
function of wavelength A and monochromatic optical depth t v . Although the model 



SolPhys_Nutto.tex; 8/03/2008; 9:14; p. 5 



C. Nutto et al. 




400 600 800 1000 1200 1400 1600 1800 2000 
Wavelength [nm] 

Figure 1 . Derivative of opacity with respect to temperature as a function of wavelength and 
the monochromatic optical depth t„. 

atmosphere extends from the photosphere to the transition region, we restrict our- 
selves to the photosphere since we are only interested in continuum radiation. Certain 
discontinuities in the plot are noticeable. The corresponding wavelengths coincide 
with the ionization wavelengths of neutral hydrogen, namely the Balmer continuum 
(A = 364. 6nm). Just above the Balmer continuum there is another ionization edge 
visible that belongs to magnesium with A = 375.6 nm. These discontinuities can be 
explained with the ionization setting in abruptly at wavelengths just smaller than 
the ionization edges. If we assume positive temperature disturbances, then higher 
temperatures give rise to a higher ionization rate which causes a greater increase of 
the opacity. Thus, the derivative of the opacity is greater for wavelengths shorter 
than the ionization edges. The effect diminishes with smaller optical depth, since 
the influence of neutral hydrogen on the opacity vanishes for higher parts of the 
atmosphere (Vernazza, Avrett, and Loeser, 1976). The influence of neutral hydrogen 
causes the steep rise of the derivative for the deeper layers. For higher layers, the 
derivative is dominated by the influence of H _ -ions (Vernazza, Avrett, and Loeser, 
1976) and the opacity is less sensitive to temperature disturbances. This can be 
explained with the combined effect of the destruction of H~ ions due to higher 
temperatures and the higher ionization of neutral hydrogen which results in a higher 
electron density. Investigations showed that if those higher electron densities are 
not taken into account the derivatives turn out to be negative for those parts of 
the spectrum where the opacity is dominated by H _ . In conclusion, the creation of 
H - due to a higher electron density wins over the destruction of H _ due to higher 
temperatures and thus the opacity increases. There is also a change of the behavior 
of the derivative for A = 1645 nm. This can be explained by the fact that above this 
wavelength the free-free processes of H~ starts to dominate the contribution to the 
opacity. 

The results for the derivative of the opacity with respect to the density are 
shown in Figure 2. First, we can point out the discontinuities at the continuum 
edges again, where the same explanation as for Figure 1 applies. For deeper layers 
the discontinuities for the Paschen continuum (A = 820.4 nm) and the Brackett 
continuum (A = 1458.4 nm) of neutral hydrogen are more pronounced than for the 
derivative of the opacity with respect to temperature. 

Knowing the derivatives of opacity, it is now possible to calculate the darkening 
and visibility functions for the most dominant solar oscillation mode with £ = 0. 
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Figure 2. Derivative of opacity with respect to density as a function of wavelength and optical 
depth. 
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Figure 3. The relative darkening functions (A„) for selected wavelengths versus position on 
the solar disk (/i = cos 9). The marks show the sampling points on the solar disk. 



4.2. Darkening and Visibility Functions 

In Figure 3 the darkening functions (A^) (see Section 2) are shown for selected 
wavelengths. Some of the these wavelengths coincide with the wavelengths examined 
by Zhugzhda, Staude, and Bartling (1996) and a comparison 2 of the results for those 
wavelengths shows that they are similar, although there are some small differences in 
the position of the change of sign. This might be caused by a small difference in the 
derivatives of the opacity derived in this paper and the derivatives used by Zhugzhda, 
Staude, and Bartling (1996) for higher layers of the atmosphere. The wavelengths 
displayed show that the position of the change of sign on the solar disk depends 
strongly on the wavelength. For the wavelengths whose opacity is mostly dominated 
by the H~-ions, e.g. 825 nm, the position of the change of sign of the darkening 
function is almost at the center of the solar disk while the function increases towards 
the limb of the Sun. Thus, for this particular wavelength the oscillation mode £ = 



2 In Zhugzhda, Staude, and Bartling (1996) the graph for A = 275 nm and A = 1645 nm should 
be switched. 
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Figure 4. The absolute visibility function versus wavelength (solid line). In addition, the 
different contributions dB v /dlnT (long dash) , d In k v jd In T (medium dash) , and d In K v /d\np 
(short dash) originating from the Planck function and the opacity are shown. 



would be best observed at the limb while it almost vanishes at the center of the 
disk. In order to get the visibility of global oscillation modes for each wavelength, 
the darkening functions have to be integrated over the solar disk. This can be done 
by either numerical integration of the darkening function for each wavelength or by 
direct evaluation of Equation (7). We choose the latter option for the calculation of 
the monochromatic visibility function. 

The result of this integration can be seen in Figure 4, where the absolute visibility 
function (5F V /{5T /Tq)) is plotted together with its different contributions from the 
Planck function and opacity. This plot clearly indicates, that the visibility functions 
are dominated by the effects of the opacity. This is due to the assumption of adia- 
batic oscillations for which the temperature and density disturbances are in phase. 
The effect of the oscillation on the opacity is big enough in order to overcome the 
disturbances of the Planck function and for shorter and longer wavelengths positive 
disturbances actually result in a darkening of the solar flux. Considering the visibility 
function, the oscillation mode £ = has the biggest influence on the disturbance of 
the flux at wavelengths around A = 800 nm. Since the visibility functions also show 
a change of sign for two wavelengths of the considered spectrum, it is interesting to 
note that for A = 480 nm and A = 1450 nm the global oscillation of mode I = 
should not be visible. This can actually be used as a test for the assumption of 
adiabatic oscillations as will be discussed below. 



5. Conclusions 

There is an essential distinction between visibility functions for the global solar 
and stellar oscillations observed by measurements of Doppler shifts and continuum 
brightness fluctuations. In the case of Doppler shifts the visibility is defined by the 
projection of the velocity vector of the oscillation modes onto the line of sight. Thus, 
it is not complicated to calculate the visibility function for every spherical harmonic. 

The derivation of the visibility functions for brightness fluctuations due to global 
oscillations is more difficult. We have succeeded in obtaining the spectral darkening 
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and visibility functions over a broad range of the visible and infrared spectrum 
of the Sun. Two simplifications of the problem make this possible: the adiabatic 
approximation and the condition 5T/Tq = constant. 

We found that the brightness fluctuations are not necessarily proportional to 
temperature fluctuations. The darkening function shows a change of sign and the 
visibility function can be negative, e.g. the brightness can decrease with the increase 
of the temperature due to radiative flux blocking. 

But there are some problems considering the visibility function. In accordance 
with Figure 4, the brightness fluctuations should be absent for 480 nm and 1450 nm. 
But this is not the case. Brightness fluctuations are observed at those wavelengths 
by SOHO/VIRGO and CORONAS/DIFOS. This clearly indicates that one of the 
assumptions is not correct, if not both of them. The condition with constant tem- 
perature disturbances should be relaxed in general. It is necessary to know the 
temperature fluctuations in dependence of the depth in the solar atmosphere. 

The second assumption in use is the adiabatic approximation. Our calculations 
confirm one of the results of Toutain and Gouttebroze (1993). If adiabatic conditions 
are assumed, the opacity perturbations show a great influence on the visibility func- 
tions. But as mentioned above, in this case the brightness fluctuations should not 
be visible for 480 nm and 1450 nm which is not in accordance with observations. 
Thus, this approximation is violated in the solar photosphere. Zhugzhda (2006) 
discovered that there are phase shifts between brightness fluctuations in different 
spectral channels of the photometric instrument DIFOS. This is only possible if 
oscillations are nonadiabatic since in this case there is a phase shift between the 
temperature and density fluctuations. Besides, the nonadiabatic p modes are cou- 
pled with temperature waves (Zhugzhda, 1983). Consequently, the darkening and 
visibility functions for nonadiabatic oscillations are complex functions. Multichannel 
observations of CORONAS together with SOHO observations make it possible to 
explore these nonadiabatic oscillations. 

The visibility and darkening functions are needed for the interpretation of the ob- 
servations of solar and stellar oscillations. The main result of the current exploration 
is that the effect of the opacity fluctuations on brightness oscillations dominates. But 
it is a rather complicated problem to find out the correct functions for adiabatic and 
nonadiabatic oscillations. The current paper just presents the first step on this road. 
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